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ABSTRACT 

We present an analysis of the evolution of 8625 Luminous Red Galaxies between 
z = 0.4 and z = 0.8 in the 2dF and SDSS LRG and QSO (2SLAQ) survey. The 
LRGs are split into redshift bins and the evolution of both the luminosity and stellar 
mass function with redshift is considered and compared to the assumptions of a passive 
evolution scenario. We draw attention to several sources of systematic error that could 
bias the evolutionary predictions made in this paper. While the inferred evolution is 
found to be relatively unaffected by the exact choice of spectral evolution model used 
to compute K+c corrections, we conclude that photometric errors could be a source 
of significant bias in colour-selected samples such as this, in particular when using 
parametric maximum likelihood based estimators. We find that the evolution of the 
most massive LRGs is consistent with the assumptions of passive evolution and that 
the stellar mass assembly of the LRGs is largely complete by z ~ 0.8. Our findings 
suggest that massive galaxies with stellar masses above 10^-^ Mq must have undergone 
merging and star formation processes at a very early stage {z > 1). This supports the 
emerging picture of downsizing in both the star formation as well as the mass assembly 
of early type galaxies. Given that our spectroscopic sample covers an unprecedentedly 
large volume and probes the most massive end of the galaxy mass function, we find 
that these observational results present a significant challenge for many current models 
of galaxy formation. 
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1 INTRODUCTION 

Luminous Red Galaxies (LRGs) are arguably some of the 
brightest galaxies in our Universe allowing us to study 
their evolution out to much higher redshifts than is possi- 
ble with samples of more typical galaxies. LRGs are known 
to form a spectroscopically homogenous population that 
can be reliably identi fied photometrically by means of sim- 
ple c olour selections (|Eisenstein et al.] l200ll : [c 



et al.l 



l2006l ). The presence of a strong 4000 Abreak in these 
galaxies also means that accurate photometric redshifts can 
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be derived for l arge samples where spect r oscopy is diffi- 
cult t o obtain (iPadmanabhan et al.l l2005l : ICoUister et al.l 
I2OO7I : lAbdalla et al.l I2OO8I ). Consequently, samples of lu- 
minous red galaxies have emerged as an importan t data 
set both for stu dies of cosmology (|Blake et all |2007| . l2008l : 
ICabre fc Gaztan aga 20 0^) and galaxy fo rmation and evolu- 
tion (|Wake et al.ll200a ICool et aLlbOOSl ). 



The study of massive galaxies in our Universe is particu- 
larly interesting since these galaxies present a long-standing 
problem for models of galaxy formation. While in the stan- 
dard ACDM cosmology, massive galaxies are thought to 
be built up through successive mergers of smaller systems 
(|White fc Reeslll978l ). many observations now suggest that 
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llGlazebrook et al. 2004; 


Cimatti et alj 20061: 


Scarlata et al.l 


l2007l:lFerreras et al.ll2009 


:|Pozzetti et al.ll200£ 


) . At the same 



time, star formation and merging activity has also been 
seen in such sys t ems between z ~ 1 and the present day 
dLin et al l l2004l: [Stanford etall l2004h and the studies of 
iBell et alJ (|2004l ') and iFaber et alJ (|200if ) have found evi- 
dence that the luminosity density of massive red galaxies has 
remained roughly constant since z~l implying a build up in 
the number density of these objects through mergers. Most 
observational results now support downsizing in star forma- 
tion - i.e most massive galaxies having lower specific star 
formation rates than their less massive counterparts - and 
many galaxy formation models are able to reproduce these 
observations e.g. by invoking quenching mechanisms suc h 



as AGN feedback (jCroton et ahlboodlPe Lucia et al.ll2006h . 
However, the more recently observed trend of downsizing in 
the mass assembly (|Cimatti et al.|[2006l : lPozzetti et al.ll2009l ) 
presents more of a challenge for galaxy formation models 
which typically predict that the most massive galaxies were 
assembled later than their less massive counterparts. 

Many of the contradictions in observational studies of 
massive galaxy formation arise principally for two reasons. 
Firstly, the way in which early-type galaxies are selected in 
different surveys can be considerably different. Morpholog- 
ical selection compared to colour selection of such objects 
will almost certainly result in different samples being cho- 
sen, particularly at high redshifts. The colour selection is 
also usually different for different samples of massive galax- 
ies as we will point out later in this paper. Secondly, many 
of the studies of massive galaxy formation mentioned so far 
are deep and narrow spectroscopic samples that will suffer 
from biases due to the effects of cosmic variance. In addition, 
one has to be careful in interpreting observational results 
as in current galaxy formation models, star formation and 
mass assembly in galaxies are not necessarily concomitant. 
So while the stars in early-type galaxies may have formed 
very early in the Universe's history, they may have formed 
in relatively small units and only merged at lower redshifts 
to create the massive ETGs we see today. 

Studies of the evolution of luminous red galaxies such 
as those analysed in this paper have already supported the 
evidence that massive galaxies have bee n passively fading 
and show very li ttle recent star formation (|Wake et al.ll200tt 
ICool et al.|[2"008l ) . In this paper, we extend this work by also 
considering the mass assembly of these systems. The 2dF 
and SDSS LRG and QSO (2SLAQ) survey presents a com- 
prehensive improvement in volume for massive galaxy sam- 
ples. The survey covers an area of 186 deg^ and reaches out 
to a redshift of 0.8 making it a promising data set to study 
the evolution of massive red galaxies. Furthermore, much of 
the 2SLAQ area is now covered by the UKIDSS Large Area 
Survey (LAS) that provides complementary data in the near 
infra-red bands for the optically selected LRGs. The advan- 
tage of near infra-red data is that the mass-to-light ratios 
and k-corrections are largely insensitive to the galaxy or stel- 
lar type and therefore the total infra-red flux in for example 
the K-band provides a good estimate of the total stellar 
mass of the galaxies. This stellar mass estimate allows us to 
study not only the star formation history but also the mass 
assembly history of these systems. As is the case for opti- 
cally selected spectroscopic samples of massive galaxies, the 



K-band selected surveys have so far bee n restricted to rela- 
tively small and deep patches of the sky (|Mignoli et al ] |2005l : 
IConselice et al.ll2007h . Clearly a large spectroscopic survey 
of massive galaxies with optical and near infra-red photome- 
try, will allow for better constraints on the evolution of these 
systems. 

In this paper, we utilise a spectroscopic sample of 
colour-selected massive red galaxies from the 2SLAQ sur- 
vey between redshifts 0.4 and 0.8. Optical photometry is 
obtained from the Sloan Digital Sky Survey supplemented 
with near infra-red data from the UKIDSS Large Area Sur- 
vey (LAS). We consider the evolution of these galaxies in 
terms of their observed colours, luminosities and comov- 
ing number densities focussing particula rly on the effects o f 
colour selection on the inferred evolution. IWake et al.1 (|2006l ) 
have presented a comprehensive analysis of the luminosity 
function of Luminous Red Galaxies using data from both 
the Sloan Digital Sky Survey and the 2SLAQ Survey. These 
authors use a subset of data from the 2SLAQ survey in the 
redshift range 0.5 to 0.6 and by comparing this galaxy popu- 
lation to a lower redshift population at 0.17 < z < 0.24 from 
SDSS, they are able to establish that the LRG LF does not 
evolve beyond that expected from a sim ple passive evolu- 
tion model at these redshifts. Meanwhile, ICool et "all (|2008l ) 
have compared the low-redshift SDSS LRG population to 
their own high-redshift sample at redshifts of ~ 0.9 and 
once again find little evidence for evolution beyond the pas- 
sive fading of the stellar populations. We extend the work of 
these authors by considering now most of the galaxies avail- 
able in the 2SLAQ data set. This enables the redshift range 
of 2SLAQ galaxies used to be broadened to 0.4 ^ z < 0.8. 
The 2SLAQ luminosity function i s presented for 862 5 LRGs 
as compared to the 1725 used bv lWake et al.l (|2006l ). These 
results are useful in filling the gap between redshift 0.4 and 
0.5 and redshift .6 and 0.8 in the IWake et al.1 (|2006l ') and 
ICool et al.l (|2008l ) data sets. The sample is split into four 
redshift bins and the evolution of the luminosity and colour 
of LRGs between these redshift bins is considered. In ad- 
dition, we also consider the stellar mass function and its 
evolution with redshift. Furthermore, we use various differ- 
ent stellar population synthesis models that are commonly 
used in the li terature, to model the LRGs, including the 
new models of iMaraston et al] (|2009l ) and quantify the sen- 
sitivity of the luminosity function estimate to changes in 
these models. The optical colours of LRGs have long been 
difficult to model using standard spectral evolution models 
jEisenstein et al.] 200ll) and this pr oblem has only recently 
been solved ( Maraston et al.ll2009l ) thereby allowing us to 
utilise the most accurate spectral evolution models of LRGs 
to date in order to infer their evolution. 

The paper is structured as follows. In §[2]we describe the 
spectroscopic 2SLAQ data set as well as SDSS and UKIDSS 
photometry for these galaxies. §[3]describes the spectral evo- 
lution models used in this paper to model the LRGs. We con- 
sider the optical luminosity function of the 2SLAQ LRGs in 
§ |4] and its evolution with redshift as well as its sensitivity 
to cosmic variance, photometric errors and changes in the 
spectral evolution models. In §[5]we present estimates for the 
LRG mass function in different redshift bins and analyse its 
sensitivity to the choice of spectral evolution model as well 
as the IMF. Finally, we discuss our results in § [6] in terms 
of models of galaxy formation and evolution. Throughout 
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this paper we assume a cosmological model with n„i=0.3, 
r2A=0.7 and h—0.7. All magnitudes are in the AB system 
unless otherwise stated. 



2 DATA 

Our dataset is a spectroscopic sample of Luminous Red 
Galaxie s from the 2dF and SPSS LRG and QSO (2SLAQ) 
survey (|Cannon et alj l2006h . This survey was conducted 
on the 2-degree Field (2dF) spectrograph on the 3.9m 
Anglo- Australian Telescope. The survey recorded spectra for 
~10000 LRGs with a median redshift of 0.55 and ~10000 
faint 2 < 3 QSOs selected from SDSS imaging data. The 
survey covers an area of 186 deg^ and extends to a redshift 
of ~ 0.8 for the LRGs. In this section, we provide a descrip- 
tion of the 2SLAQ data as well as optical and near infra-red 
photometry for these galaxies obt ained using the SDSS and 
UKIDSS LAS (jSmith et al.ll2009l ) respectively. 

In the 2SLAQ survey, the following colour and magni- 
tude cuts have been applied using SDSS DR4 photometry 
in order to is olate the main LRG population to target with 
spectroscopy (| Cannon et al. I [20061 ). 

C|| = 0.7(g -r) + 1.2(r - i - 0.18) > 1.6 (1) 
d±^ (r-i)-{g- r) /8.0 > 0.65 (2) 

17.5 ^ idev - ^ 19.8 (3) 

where g, r and i denote the dereddened SDSS model 
magnitudes, idev is the i-band deVaucouleurs mag- 
nitude and Ai is the extinction in the i-band. The 
deVaucouleurs magnitudes are obtained by fitting a 
pure deVaucouleurs profile to the 2D galaxy images. 
The definitions of both these types of magnitudes 
as well as other SDSS parameters can be found at 
http : //w ww. sdss . org/dr4/ algor ithms/photometry . html 
as well as lStoughton et al.1 (|2002h . 

The se cuts are for t he pr imary sample defined as Sam- 
ple 8 in [C annon et al.l (|2006l ) which has the highest com- 
pleteness and this is the sample that we use throughout this 
paper. M stars are expected to make up ~5% of the sam- 
ple but can be excluded by their low redshifts. In addition, 
as will be discussed later in the paper, photometric errors 
could scatter objects both in and out of the sample across 
the colo ur selection boundaries and lead to some contam- 
ination. iRoseboom et al.1 (|2006l ) have conducted a detailed 
study of the star formation histories of this 2SLAQ LRG 
sample and find that 80% of the objects, which represents 
the vast majority, are passive in nature as would be expected 
from early type systems. The rest are mainly emission line 
galaxies. Note also that the colour cuts applied in order to 
isolate LRGs are usually different for different surveys. For 
example the colour cuts used to iso late LRGs in the SDSS 
LRG sample (|Eisenstein et al.]|200ll ) results in a distribution 
that is redder than that of the 2SLAQ galaxies. If these dif- 
ferent LRG samples are to be compared therefore, further 
cuts must be applied to them to ensure a consistent colour 
selection as in IWake et al.l (|2006l ). 



In the case of the 2SLAQ LRGs, the effect of the d± cut 
is to select early-type galaxies at increasingly high redshifts 
whereas the cy cut eliminates late-type galaxies from the 
sample. Most of the sample has redshifts of 0.2 < z < 0.8 
with ~5% contamination from M-type stars jCannon et al.l 
I2OO6I ). In order to capture the main redshift distribution of 
this primary sample, only galaxies with 0.4 ^ z < 0.8 are 
used in this work. Note, however that most of the galaxies 
lie above a redshift of 0.45 and our conclusions remain un- 
changed if we choose the lower redshift limit of our sample 
to be 0.45 instead. Furthermore, we have also selected only 
galaxies with a redshift quality fiag of greater than 2. This 
results in a sample of 8625 LRGs and the redshift complete- 
ness for this sample is 76.1% (David Wake:private commu- 
nication) . All LRGs in our sample have secure spectroscopic 
redshifts. 

In Figure [TJ we plot the LRGs in the observed (g-r) ver- 
sus (r-i) plane in redshift bins of width 0.1. Isolating bright 
galaxies only with Midev < —22.5 and plotting them in the 
same colour-colour plane shows clearly how the red sequence 
is truncated in the lowest redshift bin due to the redshift de- 
pendant colour selection. The effect of this truncation on the 
luminosity function estimate will be described in more detail 
in §11 

In order to obtain accurate stellar masses for the LRGs, 
the 2SLAQ data is matched to near infra-red data from the 
UKIDSS Large Area Survey (LAS) DR5 which has a depth 
of A'vega=18.4. The galaxies are matched to within 1.1 arc- 
sees in position using the WFCAM science archiv<|3 and only 
galaxies with detections in the K-band are selected. Out of 
our total sample of 8625 2SLAQ LRGs, we have K-band data 
for 6476 of them which represents ~75% of the sample. The 
UKIDSS K-band magnitudes are in the Vega system and 
h ave been corrected to the AB system using the corrections 
of lHewett et~al] l|2006l ). 



2.1 Photometric Errors 

IWake et al.l (|2006l ) have conducted a detailed analysis of 
the effect of photometric errors on the 2SLAQ LRG sam- 
ple. These photometric errors could induce a systematic 
bias in the sample by scattering galaxies both in and out 
of the sample across the colour selection boundaries. For 
i < 19.3, the effect of the photometric errors is rela- 
tively insignificant and the colour-magnitude distributions 
of LRGs selected using single-epoch photometry and those 
selected using mul ti-epoch photometry are almost identical 
I Wake et al.ll2006h . Fainter than this however, the photo- 
metric errors may present a significant systematic bias to 
the inferred colour and luminosity evolution of LRGs in this 
survey. Note that the magnitude errors inferred using the 
single-epoch data are found to be systematically underes- 
timated. The effect of these photometric errors on the lu- 
minosity function estimate is considered in § 14.2.11 It turns 
out that different estimators of the luminosity function have 
different sensitivity to the photometric errors and this point 
is addressed in more detail in Appendix A. 
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Figure 1. The (g-r) versus (r-i) colours of 2SLAQ galaxies in the observed frame in four redshift bins - 0.4 sC z < 0.5 (top left), 
0.5 z < 0.6 (top right), 0.6 ^ z < 0.7 (bottom left) and 0.7 ^ z < 0.8 (bottom right). The dark crosses correspond to all galaxies in 
that redshift bin whereas the light dots show the brightest galaxies with Mi^^y < —22.5. 



3 K+E CORRECTIONS 

In order to calculate the rest-frame properties of LRGs 
and infer their evolution, the observed properties need 
to be transformed into the rest-frame by means of a 
k-correction. In addition, one can also correct for any 
evolutionary changes expected in the galaxy spectra by 
means of an e-cor r ection . Motivated by the work of 
IPadmanabhan et al ] (l2005[) we start by using a Pegase 
l|Fioc fc Rocca-Volmerang{ 1 19971 ) template to model the 
spectral evolution of the LRGs. In this template, the stars 
are for med in a single burst llGyr ago [z ~ 2.5). We as- 
sume a ISalpeteil (|l955l ) IMF, solar metallicity and no galac- 
tic winds or substellar objects. 

However, it has been noted by several authors that stel- 
la r population synthe s is mo dels such as Pegase and those 
of lBruzual fc CharlotI (|2003|) fail to reproduce the observed 



colours of LRGs (jEisenstein et al.ll200l|jMarastonll2005l ). Im- 
provements to these models have subsequently been made 
that involve changes to the input stellar libraries as well 
as including an additional metal-poor sub component to 
the stellar population in order to create a composite model 
(|Maraston et al.l |2009| ). Such models predict significantly 



bluer (g-r) colours compared to simple stellar population 
models and are better able to reproduce the observed broad- 
band colours of individual LRGs. For this reason, a Maras- 
ton model of age 12Gyr corresponding to a galaxy mass of 
lO^'^ Mq is also considered in this work. 

F inally, we consider the models of Char iot & Bruzual 
2007 (|Bruzual fc CharlotI [200I : lBruzua]|[2007l ) - CB07 here- 
after - with a range of different metallicities and star forma- 
tion histories in order to quantify the systematics associated 
to the different prescriptions, isochrones and stellar libraries 
of the population synthesis models. All the CB07 templates 
assume a formation redshift of zp = 3 and a lChabrieil (|2003l ) 
IMF. 

K and k+e corrections in the r-band derived from each 
of these models, are illustrated in Figure [5] In addition, k- 
corrections are also deri ved from stacked LRG spectra from 
the SDSS LRG survey (jEisenstein et al.ll200ll ) at redshifts 
0.2 and 0.5. These are very similar to each other suggest- 
ing that there is little or no colour evolution in the SDSS 
LRGs between redshifts of 0.2 and 0.5. As the 2SLAQ LRG 
spectra are not flux calibrated, a similar derivation of the k- 
corrections using these was not possible. Although the SDSS 
LRGs are selected using different colour cuts, it is still in- 
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teresting to compare their k-corrections to the ones derived 
using the various spectral evolution models and it can be 
seen that the empirical k-corrections best match those de- 
rived from a Maraston model. For this reason, the com- 
posite 12Gyr model from IMaraston et al.l |200y) has been 
used throughout this paper unless otherwise stated. We can 
also see from the right-hand panel of Figure [2] that the e- 
correction is very sensitive both to the star-formation history 
of the galaxy as well as the metallicity as galaxies with some 
residual star formation as well as those at high metallicity 
have more negative e-corrections. 

In Figure O we also plot the rest-frame (g-r) and (r-i) 
colours of LRGs derived using the different spectral evolu- 
tion models as a function of redshift. Once again it is seen 
that the Maraston model is able to reproduce the colours 
derived from SDSS composite spectra almost exactly un- 
like the Pegase model. However, we note that unlike in 
IMaraston et al.l (|2009l ). we find that including the evolu- 
tionary correction to the Maraston models, creates a mis- 
match in the colours of the model and the SDSS composite 
spectra. The SDSS observed colours are found to be red- 
der than the evolving Maraston model. We have already 
noted that the SDSS LRG colour distrib ution is redder 
than that of the 2SLAQ galaxies used in IMaraston et al.l 
(|2009l ) for comparison to their models so this difference is 
perhaps to be expected. It should also be noted that the 
observed colours derived from the spectra are derived from 
the fiber magnitudes of the galaxies. As the fibers only en- 
close light within a 3" aperture and it i s well known that 
early type galaxies have colour gradients (|Vader et al.|[l98^ : 
iFranx fc lUingworthl [l990l : iFerreras et"al]|2005l ). this could 
be why the colours appear redder but it is unlikely that 
this effect would be significant especially at redshifts of 0.5. 
Finally, the redshift 0.5 LRGs in the SDSS survey are the 
brightest an d therefore expected t o be the reddest galax ies in 
this survey (jBernardi et al.l l2005l : iGallazzi et al ] l200fih . AU 
these factors may help explain the discrepancy in the ob- 
served and model colours derived using the evolving Maras- 
ton model. 

Figure |4] shows the 2SLAQ LRGs in the observed (g-r) 
versus (r-i) colour-colour plane along with the main colour 
selection boundaries as well as the tracks produced in this 
plane by different stellar population synthesis models as well 
as SDSS composite spectra. It can be seen that both the 
Maraston as well as Pegase burst models only enter the d± 
selection at redshifts of ~0.45. Star forming galaxies with 
r > 2Gyr are clearly excluded at all redshifts by the cy 
selection. 



4 THE OPTICAL LUMINOSITY FUNCTION 

The luminosity function for 2SLAQ LRGs is calculated in 
this section in order to infer any evolution in the observed 
number density of these objects. The non-parametric Vmax 
estimator is used to calculate luminosity functions and a 
gaussian parametric form is fit by means of a minimisa- 
tion where appropriate. The gaussian parametric form pro- 
vides a better fit to the luminosity function of early-type 
galaxies compared to the more commonly used Schechter 
function and has three free parameters, 0* , M, and (Jg . The 
parametric form is given by: 



$(M) 



exp 



2^^ 



(4) 



The commonly used parametric estimator of 
ISandage et aP ([1973) (STY hereaft^er) is prone to bi- 
ases for this data set and therefore has not been presented 
in the main body of the paper. Details of the STY estimator 
and the biases induced on it are given in Appendix A. 

The Ymax estimator relies on measuring the maximum 
volume occupied by a galaxy given the survey selection cri- 
teria. By scaling the observed number density by this acces- 
sible volume, we account for the fact that only the brightest 
galaxies are observed at high redshifts in a flux-limited sam- 
ple. For the 2SLAQ survey, this maximum volume depends 
not only on the maximum redshift out to which the galaxy 
could be observed given the fiux limit of idev — Ai < 19.8 
but also the minimum redshift at which the galaxy would be 
included in the survey given the and cy colour selection 
criteria. The maximum volume for each galaxy is then given 
by: 



{z)S.^^{z)S.,^,iz)^dz (5) 



where Sd^{z), Scy (z) and Si^^^.{z) axe the selection func- 
tions for each galaxy due to the two primary colour cuts 
(Eq[T]and[2ll and the fiux limit (Eq[3]). Once this maximum 
volume has been derived for every galaxy in the sample, the 
luminosity function which is simply the number density of 
galaxies per unit brightness, is given by: 



$(M)dM = cJ2 



Vrr, 



(6) 



where c is the redshift completeness of the sample assumed 
to be 76.1% and the index, i runs over all galaxies in the 
absolute magnitude bin dA'I. 

The errors on the luminosity function are assumed to 
be Poissonian and are given by: 



a^{M)dM 



E 



1/2 



(7) 



The best-fit parametric form is derived using a fit 
to the Vmax data points. We marginalise over the normali- 
sation, (^t which is determined independantly from the ob- 
served number density of galaxies. The marginalised to 
be minimised is: 



(8) 

where $m refers to the parametric luminosity function which 
is assumed to be a gaussian function and $y represents the 
Vmax data points with errors oy- Once the best fit param- 
eters have been determined, the normalisation of the lumi- 
nosity function, 0* is then determined by matching the para- 
metric estimate of the luminosity function to the observed 
number density of galaxies using Eq. [O] 
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Figure 3. The (g-r) - left - and (r-i) -right - colours of LRGs as a function of redshift. 



Jdzfs^JdMf{M,zMM) 

where Ng is the number of galaxies, fs is the fraction of sky 
covered by the survey and f{M, z) represents the colour se- 
lection function obtained by summing the product of Sd^ (2) 
and (z) over all galaxies in the sample. The errors on the 
gaussian parameters M» and a are calculated using the la 
error ellipsoid obtained by performing the minimisation us- 
ing a Markov Chain Monte Carlo (MCMC) method. The 
main sources of error on come from cosmic variance and 
the covariance of with M«. In ij |4.2.2l we will show that for 
a large survey such as 2SLAQ, cosmic variance is unlikely to 
be a dominant systematic. The errors for cj>* quoted through- 
out the paper are therefore calculated from the la errors on 
M.. 

The dependance of the luminosity function on redshift 
is first considered and possible sources of systematic errors 
to these estimates examined later in this section. 



4.1 Redshift Evolution 

We examine the redshift evolution of the 2SLAQ LRG lu- 
minosity function in four redshift bins between redshift 0.4 
and 0.8. The luminosity function estimates are obtained us- 
ing the Vmax estimator described above. Due to the different 
absolute magnitude ranges in each redshift bin, no fits to 
the data points are shown as the inferred parameters would 
not be comparable. 

In the left-hand panel of Figure[S]we plot the luminosity 
function in four redshift bins after K correcting to redshift 
using the 12Gyr Maraston composite model. The number 
density of galaxies is found to increase slightly with redshift 
particularly in the brightest absolute magnitude bins. In the 
faintest bins there is a drop in the observed number density. 
As LRGs are known to be some of the brightest objects 
in the Universe and our survey performs an LRG selection, 
we would expect a drop in the number density of such sys- 
tems at faint magnitudes. Downsizing in the star formation 
of early type systems, is an effect that is already well ob- 
served. The fainter galaxies are expected to be less massive 
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Figure 4. The (g-r) versus (r-i) colours of 2SLAQ galaxies in the observed frame along with tracks produced by several different stellar 
population synthesis models as well as composite LRG spectra from SDSS. In the case of the models, the markers are located at redshift 
intervals of 0.1. In the case of the SDSS composite spectrum, the markers are located at redshift intervals of 0.01. 
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Figure 5. The 2SLAQ luminosity function at a rest-frame redshift after k-correction (left) and k+e correction (right) assuming passive 
evolution. 



and therefore bluer due to their higher specific star forma- 
tion rates and such objects are excluded from our sample by 
the colour selection criteria which only select out the bright- 
est and reddest objects. The Vmaa; method used to evalu- 
ate luminosity functions already accounts for those galaxies 
that may be missing from our sample due to the observa- 
tional limits set by the survey. It does not however account 
for photometric uncertainties in the sample which we study 
in detail in § 14.2.11 It is possible that genuine LRGs may 
be scattered out across the colour selection boundaries due 
to their large photometric errors rendering the sample in- 
complete at the faint end. In order to demonstrate that the 



downturn we see at the faint-end of our luminosity function 
estimates is not due to any such incompleteness, we plot 
in Figure E] the normalised histogram of the perpendicular 
distance of the faintest 2SLAQ LRGs with idev > 19 from 
the dx selection line. Most of these objects that appear faint 
and are likely to have large photometric errors, are seen to lie 
very close to the selection line - i.e. at a distance of ~0 . This 
distribution therefore suggests that there are also likely to 
be large numbers of faint galaxies on the other side of the dx 
selection. Any photometric errors would thus scatter more 
galaxies into the sample than are being scattered out result- 
ing in an overprediction in the number density of the faint 
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Figure 6. Normalised liistograms of tlie perpendicular distance 
of 2SLAQ galaxies with ijeV > 19 from the dx selection line. 



galaxies that are close to the flux limit of the survey. We 
therefore conclude that the faint-end downturn seen in our 
luminosity function estimates is genuine and not due to any 
incompleteness induced by the photometric errors. We have 
also seen in Figure [1] how the red sequence is truncated in 
the lowest redshift bin due to the redshift dependant colour 
selection. This is what leads to a drop in the number density 
in the lowest redshift bin relative to the other bins which is 
not due to any evolution in the underlying galaxy popula- 
tion. 

In the right-hand panel of Figure[5l we plot the luminos- 
ity functions after correcting for passive evolution assuming 
the same Maraston model. The luminosity functions change 
little on inclusion of this evolutionary correction suggest- 
ing it is small in these redshift ranges but in general the 
data points for galaxies with Midev < —22.5 in the different 
redshift bins now agree better. We are confident that incom- 
pleteness does not affect this bright galaxy sample and these 
results therefore suggest that this LRG population is evolv- 
ing consistently with the assumptions of a passive evolution 
model. 

4.2 Sources of Systematic Error 

In this section we examine potential sources of systematic 
error that could bias the luminosity function estimates pre- 
sented above. 

4-2.1 Photometric Errors 

In order to assess the bias that photometric errors will in- 
troduce into the luminosity function estimate, we cut our 
sample at i < 19.3 and calculate the luminosity function for 
this reduced sample of galaxies. Note that in § 12.11 we have 
already drawn att ention to the fact that previous studies 
jWake et al.ll2006l ) have shown that fainter than i = 19.3, 
galaxy samples selected through single and multi-epoch pho- 
tometry differ considerably. The single-epoch magnitude er- 
rors are also found to be systematically underestimated. 
The luminosity functions for the entire sample as well 



Figure 7. The 2SLAQ luminosity function calculated for the 
entire sample and galaxies with i < 19.3 for which photometric 
errors are thought to be insignificant along with gaussian fits to 
both. Only the bright-end of the luminosity function for the entire 
sample is shown in order to ensure that the gaussian fits are made 
over the same absolute magnitude range. 



as the i < 19.3 sample are plotted in Figure [7] along with 
gaussian fits to these derived using minimisation. Only 
the bright-end of the luminosity function for the entire sam- 
ple is shown as the sample with small photometric errors 
has no galaxies in the faintest absolute magnitude bins and 
we want to compare the two luminosity functions over the 
same absolute magnitude range. The i < 19.3 sample has 
2131 galaxies as opposed to the 8625 galaxies in the entire 
sample. The space densities from the i < 19.3 sample have 
therefore been renormalised by the fraction of galaxies in the 
complete sample that contribute at each absolute magnitude 
bin. The fit to the reduced sample with small photometric 
errors yields values of (jg = 0.45 ±0.02, M, = -21.89 ±0.05 
and 0, — 2.64 ±0 18 xl0~* where the normalisation, (fit for 
the i < 19.3 sample has also been corrected for the smaller 
space density of this sample. In this case however, a con- 
stant average correction factor is used rather than the abso- 
lute magnitude dependent corrections applied to the Vmax 
data points in order to ensure that the shape of the gaus- 
sian is maintained. The x^ fit to the entire sample on the 
other hand gives ag = 0.52 ± 0.01, M. = -21.70 ± 0.04 and 
0. = 5.51 ±g:l? xlO""*. 

As can be seen, the faintest galaxies in terms of abso- 
lute magnitude are removed from our sample on applying 
the i-band cut, as they typically have larger photometric 
errors. However, the two luminosity functions are remarka- 
bely similar in shape and there is evidence for a faint-end 
downturn even for the sample with small photometric er- 
rors. The space density at the bright-end is also lower for 
our sample with small photometric errors even after renor- 
malisation suggesting that photometric errors affect galaxies 
at all absolute magnitudes. The space density for the entire 
sample is found to be about a factor of two larger due to 
galaxies being scattered into the sample due to their large 
photometric errors. We have already shown that more galax- 
ies are likely to be scattered in across the colour selection 
boundaries than are scattered out (Figure (6]). These galax- 
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Figure 8. The 2SLAQ luminosity function calculated for the 
entire sample and in two different patches in the sky in order to 
assess the effects of cosmic variance. 



ies are bluer than more typical LRGs but appear red due 
to the photometric errors. Many of their spectra show ev- 
idence for the presence of strong emission lines suggesting 
th at they are star-forming and this has already been found 
bv lRoseboom et all (|2006 



cosmology. The results are illustrated in Figure |9] and the 
best-fit gaussian parameters to these luminosity functions 
summarised in Table [T] Both the Pegase burst and Maras- 
ton models have similar ages and star formation histories. 
The difference between them is the introduction of a metal- 
licity poor sub-component to the Maraston model and the 
use of improved stellar librarires in order to better match 
the observed colours of LRGs. However, it can be seen from 
Figure (left panel and right, top) that there is very little 
difference in the luminosity function estimate obtained from 
using these two models. 

The Pegase and CB07 models with r=lGyr produce 
very similar luminosity functions despite having different 
initial mass functions as the K-fe corrections and there- 
fore the luminosity function estimate is not very sensitive 
to the number of low mass stars. The main difference be- 
tween the Salpeter IMF assumed in the Pegase models and 
the Chabrier IMF assumed in the CB07 models, is the mass 
distribution for stars below IMq . Changing the star forma- 
tion history to have a little residual star formation shifts the 
luminosity function to slightly fainter absolute magnitudes 
and the number density in the brightest bins decreases. This 
is because including recent star formation or younger stars 
in the model means that the luminosity of the galaxy evolves 
faster and the galaxy fades quicker than in a passive evolu- 
tion scenario. 



4-2.2 Cosmic Variance 

One of the main drawbacks of the Vmax estimator is that 
it can suffer from biases due to cosmic variance. As mas- 
sive galaxies such as the ones in this sample are thought to 
be strongly clustered, the evolution in their number density 
will be very sensitive to large-scale structure. Traditionally, 
this has posed a problem for many small surveys but the 
2SLAQ survey should cover a significant enough volume to 
make this an unlikely source of problem for the LRGs being 
studied. In order to test this hypothesis, we split the LRG 
sample by declination and seperate galaxies with dec < —0.2 
and dec J5 —0.2. This results in two smaller samples with 
4340 and 4285 galaxies respectively. A luminosity function 
is then calculated for each of these subsamples as well as 
the total luminosity function for the 2SLAQ sample. These 
are plotted in Figure |8] No notable differences can be seen 
in the three luminosity function estimates confirming that 
cosmic variance is not a dominant systematic in this sample. 

4-2.3 K+e Corrections: Spectral Evolution Model 

In order to test whether the choice of spectral evolution 
model used to compute K-|-e corrections significantly affects 
the estimate of the luminosity function, we compute lumi- 
nosity functions for all the 2SLAQ galaxies between redshift 
0.4 and 0.8 using the 12Gyr Maraston composite model as 
well as a Pegase burst model, a Pegase model with an ex- 
ponentially declining star formation history with r = IGyr 
and a CB07 model with an exponentially declining star for- 
mation history with r = IGyr. Both Pegase models have an 
age of llGyr while the CB07 model has a formation red- 
shift of 3 corresponding to an age 11.35Gyr in our chosen 



4-2-4 K+e Corrections: Metallicity 

The colour-magnitude relation of early-type galaxies can be 
interpreted as a metallicity sequence an d variations in metal- 
Ucity dominate the slop e of the CMR (|Kodama fc Arimotd 
119971 : iBower et al.|[l993 ). The metallicity may therefore be 
an important parameter to vary in stellar population syn- 
thesis models when trying to model the spectral evolution 
of LRGs. In order to assess the importance of metallicity 
in the luminosity function estimate, we calculate luminosity 
functions after K-l-e correcting galaxies using a CB07 model 
at solar metallicity as well as one with variable metallic- 
ity. Both assume a formation redshift of 3 and an exponen- 
tially declining star formation history with a timescale, r 
of IGyr. In the variable metallicity model, the metallicities 
are derived from the absolute magnitudes using the colour- 
magnitu de relation of early- type galaxies in the Virgo Clus- 
ter from lBower et al.1 (119921 ) 

logio (^-^ = -0.11(M, + 21.0) - 0.18 (10) 

The effect of metallicity on the luminosity function es- 
timate is also illustrated in Figure |9l The second plot in the 
right panel illustrates the difference between the luminosity 
function estimates obtained using a CB07 model at solar 
metallicity and one with variable metallicity. The difference 
between the two estimates is found to be roughly of the same 
order as the errorbars on the individual luminosity function 
estimates. It is therefore concluded that the metallicity is 
not likely to be important in the evaluation of the luminos- 
ity function for the 2SLAQ LRGs studied here. The LRGs 
are all found to have metallicities between 0.6 and 1.5 Zq 
with a mean metallicity of roughly solar. 
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Figure 9. The 2SLAQ luminosity function computed using K+e corrections from different spectral evolution models with different 
prescriptions, isochrones and stellar libraries (left). Gaussian fits are also shown to the luminosity functions estimated using the Maraston 
12Gyr model, the Pegase llGyr burst model, a Pegase model with r=lGyr as well as a CB07 model at solar metallicity with r=lGyr. 
The best-fit gaussian parameters are summarised in Table [T] The right panel shows the difference between the luminosity function 
estimates as a function of the absolute magnitude when different parameters associated with the spectral evolution models are varied. 











Maraston 12Gyr 
Pegase Burst 11 Gyr 
Pegase T=lGyr 
CB07 T=lGyr 


0.44 ± 0.01 
0.41 ± 0.01 
0.41 ± 0.01 
0.41 ± 0.01 


-21.93 ± 0.03 
-22.04 ± 0.03 
-21.76 ± 0.03 
-21.77 ± 0.03 


3.73 ±^;^| XlO--* 
3.59 ztg-^g xlO"* 
3.58 ±[] J| xlO-'' 
3.98 ztg'lf xlO"-* 



Table 1. Summary of best-fit gaussian parameters to the luminosity function derived using different spectral evolution models 



4-2.5 K+e Corrections: Age 

In Figure[5]we also plot the Vmax estimate of the luminosity 
function after K-f e corrections to redshift assuming a Pe- 
gase single burst model at three different ages ranging from 
9 to llGyr. In the cosmological model assumed throughout 
this paper, these correspond to formation redshifts for the 
galaxy of between ~ 1.4 and 2.6. We find as expected that 
changing the formation redshift of the galaxy to later times 
results in the luminosity function moving to fainter absolute 
magnitudes or alternatively, a decrease in the number den- 
sity at a given absolute magnitude. As is the case for mod- 
els with some residual star formation, the introduction of 
younger stars into the models means that they fade quicker 



thereby shifting the luminosity function estimate to fainter 
absolute magnitudes. 



5 THE INFRA-RED LUMINOSITY FUNCTION 
AND STELLAR MASS FUNCTION 

So far we have considered the evolution of the 2SLAQ 
LRG optical luminosity function with redshift and quan- 
tified various sources of systematic error that could arise 
for this estimate. This analysis shows little evidence that 
the LRG population evolves beyond that expected from a 
simple passive evolution model between a redshift of 0.4 
and 0.8. However, such an analysis does not fully exploit 
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all the available information as the evolution of elliptical 
jalaxy populations is known to be a strong f u nction of mass 
De Lucia et a l. 2006*; Ci matti et alj I2OO6I : iFerreras et alJ 



l200g| : IPozzetti ct al. 20091. The availability of near infra- 



red data for a sizeable subsample of these LRGs means we 
can calculate accurate stellar masses for these objects from 
their K-band luminosities in order to see if the evolution of 
the number density of these galaxies is a strong function of 
the mass. 

The mass function analysis in this section is carried out 
for the subset of LRGs for which near infra-red photome- 
try is available from UKIDSS DR5. This reduced sample 
contains 6476 galaxies. As the Maraston models do not ex- 
tend to near infra-red wavelengths, the CB07 models have 
been used in this section to calculate the mass-t o-light ra- 
tios for the LRGs. Following IFerreras etal] (|2009l ). the stel- 
lar masses are computed by comparing the SDSS g,r,i,z 
fiber magnitudes with a set of 10 r-models with solar metal- 
licity, formation redshifts of 3 and exponential timescales 
between (i.e. corresponding to a Simple Stellar Popula- 
tion) and lOGyr. The CB07 models are used to generate the 
composite populations. The best-fit model, almost always 
corresponding to short formation timescales, is then used 
to determine the mass-to-light ratio in the K-heaid of the 
UKIDSS survey assuming a Chabrier or Salpeter IMF. This 
best-fit along with the K-hand petrosian magnitude allows 
us to determine the stellar mass of the galaxy. Note that no 
total model magnitudes are available in the UKIDSS sur- 
vey. However, the difference between the total model mag- 
nitudes in the SDSS i-band and the petrosian magnitudes 
in the same band are found to be of the order of 15%. Both 
magnitudes can be used as reasonable estimates for the total 
light from the galaxy. 

The mass functions are calculated as before using the 
Vmax method described in detail in §|4l The difference in the 
mass function when Ymax is determined using the best-fit 
CB07 models as opposed to a single Maraston model of age 
12Gyr, is illustrated in Figure [10] From this figure, it can 
be seen that the difference is minimal when using the two 
different spectral evolution models. Note, however that this 
may not remain the case if stellar masses could be calculated 
directly from the Maraston models as well as determining 
Vmax from it. 

The mass function for the 2SLAQ LRGs in two redshift 
bins - 0.4 ^ z ^ 0.55 and 0.55 ^ 2 ^ 0.8 - is illustrated in 
the left-hand panel of Figure 1111 These redshift limits have 
been chosen so as to ensure that incompleteness does not af- 
fect the high-mass end in either of these bins. On this plot, 
we also show the mass functions derived in two redshift bins 
for red seq uence galaxie s in th e COMBO- 17 survey (from 
Figure 9 of lBorch et aD (j2006l )). For illustrative purposes, 
we plot the K-band luminosity function in the two redshift 
bins in the right-hand panel of Figure [TT] The K-band abso- 
lute magnitudes are also determined using the best-fit CB07 
model for each LRG a nd after correcting t o the AB system 
using the correction of iHewett et all (|2006l ). 

Figure [TT] clearly shows that the 2SLAQ data set ex- 
tends to much higher stellar masses than the COMBO-17 
survey occupying a very different stellar mass range to the 
COMBO-17 red sequence and thereby allowing us to effec- 
tively probe the high-mass end of the mass function. These 
most massive galaxies are the ones expected to place the 
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Figure 10. The 2SLAQ mass function calculated using the Ymax 
method when Ymax is evaluated using a single Maraston model 
of age 12Gyr for all the LRGs and when Ymax is evaluated using 
the best-fit CB07 r model for each LRG. 



most stringent constraints on current models of galaxy for- 
mation. Our sample also probes a volume that is ~200 times 
bigger than that of COMBO-17 over the redshift range 
OA ^ z < 0.8. As is the case with the COMBO-17 data, 
there is little evidence for any evolution in the number den- 
sity at the high-mass end between a r edshift of 0.4 and 0.8 . 
This is consistent with the findings of iFerreras et al.l (|2009[ ) 
who find that the number density does not evolve in the 
highest stellar mass bins. Below ~ 10^^ Mq however, our 
sample is likely to be incomplete due to the redshift depen- 
dant selection criteria involved in isolating LRGs. 

The discrepancy between the 2SLAQ mass function pre- 
sented in th i s pap er and the COMBO-17 mass function of 
iBorch et al.l (|2006l ) for stellar masses greater than ~ 10^^ M© 
may arise due to the differing choic es of IMF used in the 
two papers. While this study uses the lChabrieJ 12003') IMF\ 
iBorch et all (|2006l ) use a truncated version of the .Salpcten 
(|l955h IMF to calculate stellar masses. In order to estimate 
how this would affect the mass function estimate, in Figure 
I12l we plot the 2SLAQ mass function over the entire redshift 
range of 0.4 ^ 2 < 0.8 calculated using both the Salpeter 
and Chabrier IMF. The Salpeter IMF pushes the mass func- 
tion towards higher masses and this results in a bigger dis- 
crepancy between the COMBO-17 and 2SLAQ data points. 
It is well known that the simple power-law form of the 
Salpeter IMF overestimates the number of low-mass stars 
in galaxies and therefore the stellar mass (jCappellari et all 
I2OO6I : iFerreras et al.l[2"008h and more r ecent analytical forms 
of the IMF such as those of iChabried (|2003h have fewer low- 
mass stars and therefore predict lower stellar masses. 

The increased number density of 2SLAQ galaxies seen 
at 3x 10"Mq to ~ IO^^Mq where it overlaps with COMBO- 
17 could be due to various reasons. Firstly, the COMBO-17 
area is much smaller than that of 2SLAQ and so it is possible 
that there could be some incompleteness at the high-mass 
end of the COMBO-17 red sequence mass function due to 
the effects of cosmic variance. More likely however, the dis- 
crepancy probably arises due to the complicated selection of 
LRGs which makes it difficult to obtain a proper estimate 
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Figure 11. The 2SLAQ stellar mass function in two redshift bins of - 0.4 < 2: < 0.55 a nd 0.55 < 2 < 0.8 a long with the COMBO-17 red 
sequence stellar mass functions in redshift bins 0.4 sC 2 ^ 0.6 and 0.6 ^ a ^ 0.8 from lBorch et al ] l|2006l '). The right-hand panel shows 
the 2SLAQ K-band luminosity function in the same two redshift bins. 
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Figure 12. The 2SLAQ mass function calculated using the Vmax 
method over the redshift range 0.4 ^ 2 < 0.8 when using a 
Salpeter and a Chabrier IMF. 



of the volume probed. As the LRGs are selected through 
colour cuts in addition to the flux limit and the colour selec- 
tion cuts the volume at low redshifts, the 2SLAQ galaxies 
will on average have a smaller Vmax than a flux-limited sam- 
ple such as COMBO-17. This will result in a larger inferred 
number density at any given mass. 



6 DISCUSSION 

In this paper we have presented luminosity and stellar mass 
function estimates for a sample of 8625 Luminous Red 
Galaxies with spectroscopic redshifts between 0.4 and 0.8 in 
the 2SLAQ survey. The evolution of the optical luminosity 
function between redshift 0.4 and 0.8 is found to be consis- 
tent with the assumptions of a passively evolving model with 
very little evidence for recent star formation. The number 
density of LRGs shows a downturn at faint magnitudes due 



to the effects of downsizing in the star formation whereby 
the less massive galaxies which are also less luminous, are 
bluer and therefore excluded from an LRG sample by the 
imposed colour selection criteria which select the reddest 
objects. There is also evidence that some bluer galaxies are 
scattered into the LRG sample due to the non-negligible 
photometric errors on them. 

The sensitivity of the optical luminosity function to 
changes in the spectral evolution models has also been stud- 
ied. The results of this analysis are summarised in Figure [131 
where we show the best-fit gaussian parameters that are fit 
to the Vmax data points for LRGs after assuming different 
star formation histories in the spectral evolution models. 
M* moves towards brighter absolute magnitudes if we as- 
sume a Pegase burst model instead of the Maraston model, 
and towards fainter magnitudes if we assume a model with 
some residual star formation signified by a star formation 
timescale of r = IGyr. The gaussians all have very similar 
widths although the Maraston model results in a gaussian 
that is slightly wider than that found with the rest of the 
spectral evolution models. 

The K-band luminosity function is calculated for about 
three quarters of the LRGs for which infra-red data is avail- 
able from the UKIDSS LAS DR5. The K-band luminos- 
ity function also shows little evidence for evolution in this 
redshift range. The K-band luminosities are used to derive 
mass-to-light ratios for the LRGs by fitting CB07 spectral 
evolution models to the multi-band photometry. These are 
then used to calculate mass functions in redshift bins. Using 
the best-fit CB07 model instead of the Maraston models in 
the estimate of the mass function, makes very little differ- 
ence to the mass function. We find that the most massive 
galaxies with M > 3 x IQ^^Mq are already well assembled 
at redshifts of 0.8 and their number density does not change 
much in the redshift range considered in this work. 

The 2SLAQ sample probes the most massive end of the 
mass function and can therefore be used to place stringent 
constraints on models of massive galaxy formation and evo- 
lution. In Figure [14] we plot the comoving number density 
of 2SLAQ galaxies as a function of redshift obtained by in- 
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Figure 13. Summary of best-fit gaussian parameters used to fit 
luminosity functions after assuming diff'ercnt spectral evolution 
models for k+e corrections. 



tegrating the mass functions presented in Figure 1111 Note 
that the horizontal errorbars simply represent the size of 
the redshift bin in our plots. We choose only galaxies with 
M > 3 X 10" Mo for which we are confident that the red- 
shift dependant selection criteria do not impose any incom- 
pleteness into the sample. We compare our estimates of the 
comoving number density with those from other surveys of 
massive galaxies as well as predictions from models of galaxy 
formation. All data points presented in this figure use stellar 
masses representative of a Chabrier IMF. Where appropri- 
ate, the quoted stellar masses in the literature have been 
corrected to this choice of IMF before plotting in Figure [TJ] 
In Table [2] we summarise the samples plotted in Figure [141 
We can see that the selection criteria for all these samples 
is varied. In the iFerreras et al.l (|2009l ) sample for example, 
there is no segregation by colour but rather a visual classifi- 
cation of early type galaxies. In this sample, no blue galax- 
ies are found above a stellar mass of 10"M^. [Pozzett i et al.l 
use different galaxy classification schemes to derive 
the galaxy stellar mass function and its evolution by galaxy 
type and find that the massive end (M> lO^^'^M©) is domi- 
nated by red spheroidal galaxies upto z~l. We can therefore 
infer that there should be a negligible fraction of blue mas- 
sive galaxies at the masses being considered in this work. We 
can also immediately see from Table [2] the unique position 
of 2SLAQ among these surveys due to its massive volume. 
This allows us to significantly reduce the size of the vertical 
errorbars in Figure [141 

From Fig ure [T^ we ca n see t h at our work agre e s wel l 
with that of iBorch et all (|2006l '): IConselice et al.l (120071 ) 
and IFerreras et al. and shows that the most mas- 

sive galaxies were already well assembled at redshift 0.8 
and there is little evide nce for their number density hav- 
ing changed since then. iFontana et"all (|2006h on the other 
hand find some evolution in the number density of mas- 
sive galaxies with redshift albeit mild up to z ~ 1.5. The 
differing comoving number densities for the different sam- 
ples can be attributed to the fact that these samples have 
all been selected in very different ways as summarised in 
Table [2] However, it is encouraging that the comoving num- 
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Figure 14. The comoving number density of massive galaxies as 
a function of redshift for the 2SLAQ survey - found by integrating 
the mass functions presented in Figure [121 between 3 X W^^Mq 
and IO'^^Mq along with the comoving number density from other 
surveys of massive galaxies as well a s mode ls of galaxy formation. 
KS06 refers to thelKhochfar fc Silkl | |2006|) model and D06 refers 
to the lDe Lucia et al.l l|2006h model. 



ber density for our colour selected sample of LRGs agrees 
well with that inferred from the visually classified sample of 
IFerreras etahl (|2009l '). 

Our observational results seem to m atch well predic- 
tions from the semi-analytical models of iKhochfar fc Silkl 
(|2006l ). These models follow the merging history of dark 
matter halos generated by the extended Press-Schechter for- 
malism. The ba ryonic physics is modelled a ccording to the 
prescriptions of IKhochfar fc BurkertI (|2005l ') and references 
therein. This model predicts that the number density of 
massive galaxies is almost constant upto redshifts of ~1 
as the number density of galaxies entering a certain mass 
bin from lower mass bins is counteracted by the number 
density of galaxies leaving that mass bin for a higher mass 
bin. If this were the case, we would expect the number den- 
sity of galaxies in the highest mass bin to decrease as a 
function of redshift even if it is constant in the interme- 
diate mass bins. We compute the comoving number den- 
sity in the higher mass bin of 12 < logj^Q(A/*/A/0) < 12.5 
and in our two redshift bins - 0.4 ^ z < 0.55 and 0.55 ^ 
z < 0.8. The number densities per unit dex in stellar mass 
thus obtained are log(n)//i^Mpc"^dea::"^ = -5.39±o;}3 and 
\og{n)/h^Mpc~^dex~^ = -5.55±o;o9 respectively. There- 
fore there is no evidence that the comoving number density 
of the most massive galaxies is decreasing with increasing 
redshift. The ag reement between our data points and the 
IKhochfar fc Silk ( 2006 ) models shown in Figure [14] should 
therefore be treated with caution as the constant number 
density of massive galaxies in the models, is only maintained 
due to a fine-tuning of objects between mass bins and we find 
no evidence that this fine-tuning actually occurs. 

Our results also do not match at all those gener- 
ated from models b ased on the Millennium simulations 
(|De Lucia et al.|[2006l ). In these models, AGN feedback is in- 
voked to shut off star formation after a characteristic mass 
scale in order to reproduce the observed colour-bimodality 
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Rcdshift Range 


Area 


Number of Galaxies 


Selection 


2SLAQ1 
COMBO-172 
Palomar/DEEP2=* 
GOODS-MUSIC" 
GOODS Massive Galaxies^ 


0.4 - 0.8 
0.2 - 1.0 
0.4 - 2.0 
0.0 - 4.0 
0.4 - 1.2 


186.2 deg^ 
1.0 dcg^ 
1.36 deg2 
143.2 arcmin-^ 
320 arcmin'^ 


8625 
-25,000 
4571 
2931 
457 


Luminous Red Galaxies 
All Galaxies 
K-band selected massive galaxies 
K-band selected red sequence 
Visually classified ETGs 



Table 2. Summary of M assive Galaxy Sample s f or which stellar mass functions have been calculated ■ 
^Conselice et al.l |20o3), iFontana et al.l l|2006l '). iFerreras et al.l |20oi). 



^This Work, feorch et ahl 1I2OO6I) . 



and luminosity function at redshift 0. As a consequence, the 
downsizing in star formation is reproduced but the growth of 
massive galaxies at high redshifts of greater than 0.5 is pro- 
hibited and dry mergers become the most dominant mech- 
anism for mass growth in such galaxies. This is clearly in- 
consistent with our observations which suggest not only that 
significant star formation in massive galaxies has ceased at a 
redshift of 0.8 but also that these galaxies were already well 
ass embled at these redsh ifts. The main difference be tween 
the lDe Lucia et al.1 (|2006l ') and lKhochfar fc Silkl (120061 ) mod- 
els is the feedback with respect to AGN activity rather than 
the fact that the former follows a numerical simulation for 
the evolution of the dark matter halos whereas the latter 
relies on an analytical approach. If cooling is suppressed, 
the efficiency of massive galaxy formation at early epochs 
is much lower accounting for the lower number densities in 
the De Lucia et al.l (l2006h models. We also note that the 

Lucia et al. I (|2oodr models are computed from the Mil- 
lennium sim ulations, which use a c osmology with (78=0.9, 
whereas the iKhochfar fc Siu3 ()2006h models have been up- 
dated fo r this work with the later WMAP5 cosmology with 
(78=0.8 l|Dunklev et aLlbOOgl ). As the AGN activity affects 
the star formation history of the galaxy, we conclude that 
the star formation and quenching mechanisms invoked in 
the different galaxy formation models need to be revisited 
in light of our new constraints on the very high mass end of 
the stellar mass function. 

iGranato et al.l (|2004l ) on the other hand have proposed 
a way in which to bring about the anti-hierarchical forma- 
tion of the baryonic component of galaxies in models while 
still working within the framework of the ACDM cosmol- 
ogy. This is essentially done by invoking feedback from su- 
pernovae as well as the nuclear activity in massive galax- 
ies. These processes slow down star formation in the least 
massive halos and drive gas outfiows thereby increasing the 
stellar to dark matter ratio in the more massive halos and 
ensuring that the physical processes acting on baryons are 
able to effectively reverse the order of formation of galaxies 
compared to dark matter halos. These models have already 
been shown to matc h the local K-band l uminosity function 
of massive galaxies (|Granato et al. as well as the lu- 

minosity functio n at redshifts of ~ 1.5 ( Granato et al.|[2004l : 
ISilva et al. I l2005h . Reproducing our observational results for 
the most massive galaxies at intermediate redshifts would 
therefore be an important test for these models. 

Recently iDekel et al.l (|2009l ) have also proposed an al- 
ternative galaxy formation model whereby bursts of star for- 
mation can occur in massive galaxies without the need for 
violent merging events. These authors suggest that the mas- 
sive galaxies, which reside in the centers of filaments, are fed 
by cold gas streams that penetrate the shock heated media 



of massive dark matter halos thereby inducing star forma- 
tion. 

Given the small errorbars in our estimates of the co- 
moving number density compared to previous studies and 
the fact that our sample has secure spectroscopic redshifts 
for all objects and a massive area that reduces errors due 
to cosmic variance, this new observational result presents a 
significant challenge for many current models of galaxy for- 
mation and requires a reinvestigation of the feedback mech- 
anisms involved. 



7 CONCLUSION 

This paper has examined the evolution of the optical and 
near infra-red luminosity function as well as the stellar mass 
function of 8625 LRGs from the 2SLAQ survey between red- 
shift 0.4 and 0.8. We have demonstrated the unique position 
of 2SLAQ LRGs among other massive galajcy samples due to 
the large volume probed by the sample as well the availabil- 
ity of spectroscopic redshifts for all galaxies and near infra- 
red data for ~75% of them. This has allowed us to probe 
the very high-mass end of the stellar mass function which 
places the most stringent constraints on models of galaxy 
formation . We hav e also used the spectral evolution models 
of Maras ton et al.l (|2009i ) which are the most accurate mod- 
els of LRGs to date to study the evolution of these objects 
although our conclusions change little if other models are 
used instead. Specifically we draw the following conclusions: 

• The evolution of the optical luminosity function of 
LRGs between redshifts 0.4 and 0.8 is consistent with the 
assumptions of a passive evolution model where the stars 
were formed at high redshift with little or no evidence for 
recent episodes of star formation. 

• A downturn is seen at the faint end of the LRG lu- 
minosity function due to the effects of downsizing in the 
star formation. The faintest galaxies are expected to be less 
massive and therefore bluer than the bright objects and such 
objects are excluded from the LRG sample due to the im- 
posed colour selection of these objects which selects only the 
reddest galaxies. 

• We find that photometric errors may induce a signifi- 
cant bias into our sample of LRGs and scatter galaxies both 
in and out of our sample across the colour selection bound- 
aries. We show that more galaxies are likely to be scattered 
into the sample than are scattered out leading to an over- 
prediction of the observed space density but the shape of 
the luminosity function does not change if we remove the 
galaxies with large photometric errors. 
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• The LRG luminosity function is found to be unaffected 
by cosmic variance due to the large volume occupied by the 
sample. 

• The luminosity function is also relatively insensitive to 
the choice of spectral evolution model and the metallicity 
although models with some residual star formation shift the 
luminosity function towards fainter absolute magnitudes. 

• The stellar mass function for these LRGs for M > 
3 X 10^^ Mq also shows little evidence for evolution between 
redshifts 0.4 and 0.8 suggesting that these most massive sys- 
tems were already well assembled at redshifts of 0.8. This 
is consistent with the emerging picture of downsizing in the 
mass assembly of massive galaxies. 

• The stellar mass function estimate is also relatively in- 
sensitive to the choice of spectral evolution model assumed 
in the calculation of Vmax- However, different choices of the 
stellar initial mass function will shift the mass function es- 
timate as found in previous studies. 

• The comoving number density of LRGs with AI > 
3 X 10^^ Mq has changed little between redshifts 0.8 and 0.4. 
The same is true for LRGs with M > 10^^ Mq. This is con- 
sistent with other observational results for massive galaxy 
samples and does not agree with the predictions of most 
curre nt galaxy form a tion models. We find that the mod- 
els of iGranato et all ()2004h may be promising in matching 
our observations although they are yet to be compared with 
massive galaxy samples at intermediate redshifts. 

Overall, our results support the emerging picture of 
downsizing in both the star formation and mass assembly 
of early type galaxies and present a significant challenge for 
current models of galaxy formation. We have shown our find- 
ings to be robust to changes in spectral evolution models, 
cosmic variance and photometric errors and conclude that 
these new observations of the most massive and most lumi- 
nous galaxies in our Universe will need to be reconciled with 
the models if progress is to be made in the field of massive 
galaxy formation and evolution. 
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APPENDIX A: THE STY ESTIMATOR 

The most commo n parametric e s timat or of the luminosity 
function is that of lSandage et al.] (| 19791 ) (STY). A paramet- 
ric form is assumed for the luminosity function e.g. a gaus- 
sian as in this paper. The parameters, a are then solved for 
by maximising the likelihood: 



^(a) = nP' 



(Al) 



where pi is the probability of a galaxy with absolute magni- 
tude, M being selected in the flux limited, colour cut sample 
and is given by Eq[ 
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Figure Al. The non-parametric Vmax data points for the LRGs 
in the entire sample along with the STY gaussian fit to the entire 
sample. 
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where f{Mi,Zi) is the apparent magnitude incompleteness 
function for objects with absolute magnitude. Mi and red- 
shift, Zi. This incompleteness occurs due to various reasons. 
Firstly, not all galaxies target ted in a survey wi ll have red- 
shifts determined succesfuUy. IWake et al.l (120061 ) have con- 
sidered this incompleteness function for 2SLAQ LRGs using 
carefully constructed 2SLAQ survey masks. Although they 
find a slight dependance of this incompleteness function on 
the (g-r) colour, including this functional dependance as op- 
posed to considering a constant incompleteness of 76.1% in- 
troduces less than a 1% change in their luminosity function 
estimate. We therefore assume a constant redshift incom- 
pleteness of 76.1% that is independant of both the redshift 
and the absolute magnitude. 

In the case of a colour-selected sample such as 2SLAQ 
however, there may be additional sources of incompleteness 
as the space density of these galaxies is a function of both 
the flux limit and the colours. The 2SLAQ LRGs at different 
redshifts occupy different distributions in the colour absolute 
magnitude plane due to the redshift dependant colour selec- 
tion. This colour selection could introduce an incomplete- 
ness function that is both a function of absolute magnitude 
as well as redshift and this selection function has to be ex- 
plicitly included in the integral in Eg IA2I We have checked 
that using the absolute magnitude averaged selection func- 
tion rather than the 2d function f{Mi, Zi), makes very little 
difference to the STY estimate of the luminosity function. 
In this case, the galaxy incompleteness can be taken outside 
the integral and the likelihood i n the STY method a djusted 
directly e.g. (|Zucca et al.ll 19941 : IWillmer et al.ll2006D . 



(A3) 



where the weights associated with each individual galaxy 
are Wi — 1/S{mi). We initially calculate luminosity func- 



tions assuming Wi = 1 for the entire sample and this is 
shown by the dashed hne in Figure lAll This clearly illus- 
trates that the Ymax and STY estimators differ considerably 
from each other at the faint-end. The STY fit parameters are 
M. = -21.32, a = 0.65, (p, = 9.97 x 10~* for the entire sam- 
ple and these values are very different to those derived from 
the x^"fits in § m It has already been noted in this paper 
that fainter than i ~ 19.3, photometric errors may become 
significant in scattering LRGs both in and out of the colour 
selection boundaries. This analysis therefore suggests that 
the STY estimator is particularly sensitive to photometric 
errors for colour selected samples and the weights Wi start 
to differ significantly from 1 for galaxies with large photo- 
metric errors. In order to calculate these weights, one would 
have to do a Monte Carlo simulation of galaxies assuming a 
certain spectral evolution model, add noise to this simula- 
tion and then look at the completeness function Sirrii) for 
each galaxy by considering which galaxies are sel ected on ap- 
plying the 2SLAQ colo ur selection criteria e.g. (| Fried et al.l 
I2OOII : (Wolf et all 120031 ). This is the subject of future work. 
The STY estimator in its traditional form is clearly biased 
for a colour-selected sample with non-negligible photometric 
errors and for this reason, it has not been presented in the 
main body of this paper. We note that when we move to 
high redshift bins where the galaxies are generally brighter 
and the photometric errors smaller, the STY estimator and 
Vmax estimator agree very well as it is clearly reasonable in 
these regimes to assume Wi — 1. 
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